Implicit Time Integration Scheme Research Articles

A phase field-immersed boundary-lattice Boltzmann coupling method for fluid-structure interaction analysis

Fluid-structure interaction (FSI) is critical in numerous engineering and scientific applications. This paper presents a numerical model for FSI simulations that integrates the lattice Boltzmann method (LBM) for fluid dynamics, the phase field method (PFM) for structural deformation and fracture, and the immersed boundary method (IBM) for fluid-solid coupling. The fluid employs an explicit time integration scheme with a smaller time step than the solid, which uses an implicit time integration scheme, satisfying the CFL condition. During the synchronization of fluid-solid interaction information, a multi-substep strategy is applied. The IBM is implemented using the iterative velocity method for boundary treatment. This approach facilitates the information transfer between domains with different time steps, enhancing the coupling accuracy through iterative computation. The effectiveness of the improved immersed boundary-lattice Boltzmann (IB-LBM) has been verified by the classical cylindrical flow case, which effectively prevents the particle penetration phenomenon and improves the calculation accuracy. The phase field-immersed boundary-lattice Boltzmann method (PF-IB-LBM) is subsequently deployed to simulate the deformation of an aluminum plate under hydrostatic pressure and the aftermath of a dam-break impact on an elastic plate. Verification of the coupling method substantiates its feasibility and efficacy. In addition, the proposed PF-IB-LBM is implemented to simulate the entire process of structural deformation and failure in response to fluid dynamics, which further confirms its capability to predict discontinuous mechanical behaviors such as damage and crack propagation under intricate fluid flows. The results emphasize the flexibility and efficiency of the proposed approach, providing a robust framework for solving more complex scenarios in engineering FSI problems.

Implicit direct time integration of the heat conduction problem in the Method of Matched Sections

The paper is devoted to further elaboration of the Method of Matched Sections as a new branch of finite element method in application to the transient 2D temperature problem. The main distinction of MMS from conventional FEM consist in that the conjugation is provided between the adjacent sections rather than in the nodes of the elements. Important feature is that method is based on approximate strong form solution of the governing differential equations called here as the Connection equations. It is assumed that for each small rectangular element the 2D problem can be considered as the combination of two 1D problems – one is x-dependent, and another is y-dependent. Each problem is characterized by two functions – the temperature, , and heat flux . In practical realization for rectangular finite elements the method is reduced to determination of eight unknowns for each element – two unknowns on each side, which are related by the Connection equations, and requirement of the temperature continuity at the center of element. Another salient feature of the paper is an implementation of the original implicit time integration scheme, where the time step became the parameter of shape function within the element, i.e. it determines the behavior of the Connection equations. This method was early proposed by first author for number of 1D problem, and here in first time it is applied for 2D problems. The number of tests for rectangular plate exhibits the remarkable properties of this “embedded” time integration scheme with respect to stability, accuracy, and absence of any restrictions as to increasing of the time step.

Adjacency-based, non-intrusive model reduction for vortex-induced vibrations

Vortex-induced Vibrations (VIV) pose computationally expensive problems of high practical interest to several engineering fields. In this work we develop a non-intrusive, reduced-order modeling methodology, applied to two-dimensional (2D), VIV cases subject to a laminar, incompressible flow. Performing a physics-informed approximation of the Arbitrary Lagrangian–Eulerian (ALE) incompressible Navier–Stokes (NS), we derive a discrete-time, quadratic-bilinear model structure for the velocity flowfield on a reference domain. This structure, along with a predefined sparsity pattern motivated by the adjacency-based sparsity of the discretized NS-ALE operators, leads to the formulation of a sparse, full-order model (sFOM) inference problem. Thus, the data-driven inference task requires solving many “local” least squares (LS) problems, isolating the contribution of geometrically “nearest neighbours” for each degree of freedom. Numerical aspects of inference such as data centering and regularization, as well as the direct enforcement of boundary conditions under the sFOM formulation are extensively discussed. The inferred, sFOM operators are then projected to a non-intrusive reduced-order model (ROM) for the velocity flowfield via the Proper Orthogonal Decomposition (POD). The resulting non-intrusive ROM (sFOM-POD) is coupled with the first-principle, 2D solid oscillation equations, resulting to a hybrid physics-informed/first-principle VIV dynamics model, simulated using an implicit time integration scheme. The mapping of the coupled solution from the ALE reference domain to the current configuration is also presented. This methodology is applied to two testcases of an elliptical, non-deformable solid mounted on springs, subject to Re=90,180 flows. Numerical results indicate a successful coupling between the data-driven flowfield and solid dynamics, with prediction errors of less than 3% for both the flowfield and the solid oscillation. A comparative study with respect to the sFOM-POD dimension indicates the robustness and potential of the approach.

A highly efficient computational approach for fast scan-resolved simulations of metal additive manufacturing processes on the scale of real parts

This article proposes a novel high-performance computing approach for the prediction of the temperature field in powder bed fusion (PBF) additive manufacturing (AM) processes. In contrast to many existing approaches to part-scale simulations, the underlying computational model consistently resolves physical scan tracks without additional heat source scaling, agglomeration strategies or any other heuristic modeling assumptions. A growing, adaptively refined mesh accurately captures all details of the laser beam motion. Critically, the fine spatial resolution required for resolved scan tracks in combination with the high scan velocities underlying these processes mandates the use of comparatively small time steps to resolve the underlying physics. Explicit time integration schemes are well-suited for this setting, while unconditionally stable implicit time integration schemes are employed for the interlayer cool down phase governed by significantly larger time scales. These two schemes are combined and implemented in an efficient fast operator evaluation framework providing significant performance gains and optimization opportunities. The capabilities of the novel framework are demonstrated through realistic AM examples on the centimeter scale including the first scan-resolved simulation of the entire NIST AM Benchmark cantilever specimen, with a computation time of less than one day. Apart from physical insights gained through these simulation examples, also numerical aspects are thoroughly studied on basis of weak and strong parallel scaling tests. As potential applications, the proposed thermal PBF simulation approach can serve as a basis for microstructure and thermo-mechanical predictions on the part-scale, but also to assess the influence of scan pattern and part geometry on melt pool shape and temperature, which are important indicators for well-known process instabilities.

A numerical framework for simulating progressive failure in composite laminates under high-cycle fatigue loading

In this work, a recently proposed high-cycle fatigue cohesive zone model, which covers crack initiation and propagation with limited input parameters, is embedded in a robust and efficient numerical framework for simulating progressive failure in composite laminates under fatigue loading. The fatigue cohesive zone model is enhanced with an implicit time integration scheme of the fatigue damage variable which allows for larger cycle increments and more efficient analyses. The method is combined with an adaptive strategy for determining the cycle increment based on global convergence rates. Moreover, a consistent material tangent stiffness matrix has been derived by fully linearizing the underlying mixed-mode quasi-static model and the fatigue damage update. The enhanced fatigue cohesive zone model is used to describe matrix cracking and delamination in laminates. In order to allow for matrix cracks to initiate at arbitrary locations and to avoid complex and costly mesh generation, the phantom node version of the eXtended finite element method (XFEM) is employed. For the insertion of new crack segments, an XFEM fatigue crack insertion criterion is presented, which is consistent with the fatigue cohesive zone formulation. It is shown with numerical examples that the improved fatigue damage update enhances the accuracy, efficiency and robustness of the numerical simulations significantly. The numerical framework is applied to the simulation of progressive fatigue failure in an open-hole [±45]-laminate. It is demonstrated that the numerical model is capable of accurately and efficiently simulating the complete failure process from distributed damage to localized failure.

A comparative study of fully implicit staggered and monolithic solution methods. Part II: Coupled excitation–contraction equations of cardiac electromechanics

A comparative study of fully implicit staggered and monolithic solution algorithms for the coupled bidomain equations of cardiac electrophysiology was discussed in Part I (JCAM 407:114021, 2022). In this follow-up study, we further employ the similar approach and present an extensive comparative analysis between the staggered and monolithic solution techniques for the coupled partial differential equations (PDEs) of cardiac electromechanics both in the monodomain and bidomain setting. For the monolithic solution technique, we adopt the work of Göktepe et al. (CM 45: 227-243, 2010) and Dal et al. (CMAME 253: 353-336, 2013), where the governing PDEs of the excitation–contraction problem are solved simultaneously and all rate equations are treated with an implicit time integration scheme. For the staggered scheme, however, we simply suggest a decoupled solution of the governing PDEs while keeping other aspects of the monolithic algorithm the same, such as utilization of the implicit time integration scheme and coupling of ordinary differential equations, which describe the state variables of cardiac electrophysiology, to the parabolic PDE of the electrophysiology. Both solution algorithms are applied to several problems, where we simulate regular planar wave propagations and scroll waves for different spatial and temporal resolution and material parameters. The comparison between the solution schemes is performed in terms of accuracy, efficiency and stability. We reveal that for regular wave propagation the results of different solution algorithms are identical. The only deviation is observed during scroll wave propagation, if the parameters controlling the mutual interaction between the mechanical and electrical field are relatively large and at the same time if the time increment is big, e.g., Δt=2.0ms. Nevertheless, the suggested staggered solution scheme yields almost identical results even for high coupling if smaller time increments are used. Besides, the staggered scheme provides an enormous gain in computational efficiency, particularly, as the number of degrees of freedom is increased and we do not encounter any stability issue arising specifically from the decoupled solution of the governing PDEs. Similar to the conclusion that was made for cardiac electrophysiology in Part I, we report that the suggested fully implicit staggered solution algorithm has a tremendous application potential for computer simulations of excitation–contraction coupling of the heart both in the monodomain and bidomain setting.

A hybridizable discontinuous Galerkin formulation for the Euler–Maxwell plasma model

This work introduces a hybridizable discontinuous Galerkin formulation for the simulation of ideal plasmas governed by the Euler–Maxwell equations. The approach is based on a monolithic source-based coupling of the fluid and electromagnetic subproblems, along with a fully implicit time integration scheme, and a projection-based divergence correction method to enforce the Gauss laws in challenging scenarios. The numerical examples demonstrate the high-order accuracy, efficiency, and robustness of the proposed formulation and validate it against problems of increasing complexity, ranging from single-physics cases to weakly and fully coupled electromagnetic plasma simulations.

High-order composite implicit time integration schemes based on rational approximations for elastodynamics

A new approach for developing implicit composite time integration schemes is established starting with rational approximations of the matrix exponential in the solution of the equations of motion. The rational approximations are designed to have the same effective stiffness matrix in all sub-steps. An efficient algorithm is devised so that the implicit equation for a sub-step is in the same form as that of the trapezoidal method. The proposed M-schemes are Mth order accurate using M sub-steps. The amount of numerical dissipation is controlled by a user-specified parameter ρ∞. The (M+1)-schemes that are (M+1)th order accurate using M sub-steps can also be constructed, but the amount of numerical dissipation is built in the schemes and not adjustable. The order of accuracy of all the schemes is not affected by external forces and physical damping. Numerical examples demonstrate that the proposed high-order composite schemes are effective for introducing numerical dissipation and allow the use of large time step sizes in wave propagation problems. The source code written in MATLAB is available for download at: https://github.com/ChongminSong/CompsiteTimeIntegration.

Development of a Composite Implicit Time Integration Scheme for Three-Dimensional Discontinuous Deformation Analysis

Discontinuous deformation analysis (DDA) is a discontinuum-based and implicit method for investigating the deformational behavior of block systems. The constant acceleration integration (CAI) scheme characterized by unconditional stability is employed in the traditional DDA. In this study, the problems of the CAI scheme regarding the time step in DDA are pointed out. A too large or too small time step size adopted in the CAI scheme will have adverse effects on the DDA computation. To overcome the weaknesses, an alternative composite implicit time integration (CITI) scheme, which is a combination of the trapezoidal rule and the three-point backward Euler method, is implemented in the three-dimensional (3D) DDA method. Verification examples and slope numerical simulations are presented to illustrate the accuracy and effectiveness of the proposed methodology. The results showed that the CITI scheme can overcome the numerical error caused by the large time step size, and the algorithm damping is closely related to the time step size and the selected splitting parameter.

An adaptive dynamic phase-field method using the variable-node elements for cohesive dynamic fracture

This work presents an adaptive phase-field method incorporated into a finite element framework combined with variable-node elements to investigate the cohesive dynamic fracture. The proposed framework builds upon the regularized phase-field cohesive-zone model as its foundation utilizes a hybrid form of the history field to drive the crack evolution. A staggered iteration scheme is utilized to compute the displacement and phase-field variables for the coupled elastic-phase field system. The error indicator, utilizing the phase-field as well as the history strain variables, is employed to control the adaptive refinement process. The variable-node element technique facilitates adaptive mesh refinement, which is simple and flexible to act as a transition element between coarse and refined elements. In addition, an implicit HHT time integration scheme is employed for temporal discretization. Several standard problems are presented to showcase the efficiency and accuracy of the proposed method, highlighting its superiority when compared to the results documented in the literature. The results show that the proposed framework can significantly improve computational efficiency without affecting the accuracy of numerical results.

The edge-based smoothed FEM with [formula omitted]-Bathe implicit temporal discretization scheme for the analyses of underwater wave propagation problems

The underwater acoustic wave propagation problem is an important issue in ocean engineering field. The classical finite element approach with edge-based smoothed gradient smoothing technique is incorporated into the novel ρ∞-Bathe implicit direct time integration scheme in this paper for solving underwater transient acoustic wave propagation problems. This edge-based smoothed FEM (ES-FEM) and ρ∞-Bathe method are employed to perform the space and time domain discretization, respectively. Due to the appropriate softening effects from the ES-FEM, the space discretization errors can be markedly reduced. Meanwhile, the possible numerical error related to the temporal discretization also can be significantly suppressed by the novel ρ∞-Bathe method which possesses the controllable numerical damping effects. Through detailed analysis of the total dispersion error, we find that the numerical results from the proposed scheme can converge monotonically to the exact solutions when employing the decreasing temporal discretization intervals and the possible numerical anisotropy effects are also effectively suppressed. These excellent numerical performance clearly distinguish the present method from the conventional approaches and make it to be particularly suitable for complicated acoustic wave propagation analysis. Finally, various typical numerical experiments are performed to illustrate the capacities of the proposed scheme in tackling underwater transient acoustic wave propagation problems.

Numerical modeling of wave propagation phenomena in thermo-poroelastic media via discontinuous Galerkin methods

We present and analyze a high-order discontinuous Galerkin method for the space discretization of the wave propagation model in thermo-poroelastic media. The proposed scheme supports general polytopal grids. Stability analysis and hp-version error estimates in suitable energy norms are derived for the semi-discrete problem. The fully-discrete scheme is then obtained based on employing an implicit Newmark-β time integration scheme. A wide set of numerical simulations is reported, both for the verification of the theoretical estimates and for examples of physical interest. A comparison with the results of the poroelastic model is provided too, highlighting the differences between the predictive capabilities of the two models.

High-order implicit time integration scheme with controllable numerical dissipation based on mixed-order Padé expansions

A single-step high-order implicit time integration scheme with controllable numerical dissipation at high frequencies is presented for the transient analysis of structural dynamic problems. The amount of numerical dissipation is controlled by a user-specified value of the spectral radius ρ∞ in the high frequency limit. Using this user-specified parameter as a weight factor, a Padé expansion of the matrix exponential solution of the equation of motion is constructed by mixing the diagonal and sub-diagonal expansions. An efficient time-stepping scheme is designed where systems of equations, similar in complexity to the standard Newmark method, are solved recursively. It is shown that the proposed high-order scheme achieves high-frequency dissipation, while minimizing low-frequency dissipation and period errors. The effectiveness of the provided dissipation control and the efficiency of the scheme are demonstrated by numerical examples. A simple guideline for the choice of the controlling parameter and time step size is provided. The source codes written in MATLAB and FORTRAN are available for download at: https://github.com/ChongminSong/HighOrderTimeIntegration.

A hybrid implicit-explicit discontinuous Galerkin spectral element time domain (DG-SETD) method for computational elastodynamics

SUMMARY A hybrid implicit-explicit (IMEX) discontinuous Galerkin spectral element time domain (DG-SETD) algorithm is proposed to simulate 3D elastic wave propagation in inhomogeneous media. In this method, the original problem can be divided into a number of well designed subdomains, and the mesh generation of different subdomains is completely independent, thus allowing flexible spatial discretization of complex geometry. The neighboring subdomains are connected by a Riemann transmission condition (RTC), and spectral elements with different interpolation orders can be used in different subdomains to maximize the computational efficiency of multiscale problems to facilitate parallel computing for different subdomains. In particular, the explicit or implicit time iteration scheme can be appropriately selected for a subdomain according to the size of its mesh elements to increase the global time step increment, thus giving higher computational efficiency: For subdomains with coarse meshes, the explicit time integration scheme is used and the time step increment is limited by the Courant−Friedrichs−Lewy (CFL) stability condition; for subdomain with fine structures and thus fine meshes, an implicit time integration scheme is used so that a large time step increment can be used without affecting the stability. In addition, a jump condition of displacement and velocity is introduced to accurately simulating fractures and faults, including lossless and viscous fractures with plane, curve or cross structures. This avoids the volume modeling of the extremely thin fracture structures, and effectively reduces the number of degrees of freedoms (DoFs) in the discretized system without the loss of accuracy. The accuracy, robustness and efficiency of the DG-SETD algorithm are demonstrated by quantitative comparisons of the waveforms with the commercial finite element software COMSOL.

An efficient SBFEM-based approach for transient exterior vibro-acoustic analysis of power-law functionally graded shells

An efficient approach for transient exterior vibro-acoustic analysis of power-law functionally graded (FG) shells is developed using the scaled boundary finite element method. In the structural formulation, a shell element is treated as a three-dimensional continuum and its middle surface is represented with a quadrilateral spectral element. Along the thickness, the middle surface is scaled and the displacements are approximated using a quadratic Lagrange interpolation. The assumed natural strain method is applied to treat numerical locking and the integral along the thickness is performed analytically. In the acoustic formulation, velocity potential acts as the basic unknown. An infinite fluid is truncated by a spherical surface. The exterior field is simulated through the improved doubly asymptotic open boundary while the interior region is split into subdomains with their impedance evaluated by the improved continued-fraction method. The structure and fluid are discretized independently. The collocation procedure is utilized to perform the data transfer across the nonconforming interface. A high-order implicit time integration scheme is applied to solve the coupled system of equations. Numerical examples demonstrate that the proposed approach is stable, accurate and highly efficient. The effects of the geometric and material parameters on the vibro-acoustic behaviors of FG shells are systematically studied.

A reduced-order cutFEM approach to model complex moving sound sources

While traditional finite element based vibro-acoustic analysis of structures is done on stationary structures, in many occasions the structure under investigation is actually moving. To perform accurate finite element analysis on the vibro-acoustic performance of such structures requires the evaluation of a time-varying system, leading to several challenges. Of paramount importance are a strategy to update the mesh to the new configuration and an adequate time integration scheme that does not introduce large numerical artifacts. Additionally, the time-varying nature leads to an increase in computational complexity. In this paper a modeling strategy based on cut finite elements (cutFEM) is introduced that aims to tackle these challenges. The cutFEM approach can be seen as a fictitious domain approach in which the source is moving over a static mesh, cutting through elements where the actual source is located, using a level set description. The key advantage of this approach is that no remeshing and/or mesh morphing is required. It is shown how both implicit and explicit time integration schemes can be used to propagate the solution while the source is moving. Additionally, several model order reduction strategies are compared to speed up the online evaluation of the reduced order model.

An adaptive numerical method for simulating diffusion flame jets

In this paper, we present an adaptive numerical methodology for a reacting three-dimensional low-Mach number flow. This computational strategy combines an adaptive mesh refinement (AMR), an implicit–explicit time stepping strategy (IMEX), an extension of increment-pressure projection method, and a mixture fraction to model the chemistry combustion dynamics. To accurately resolve sharp gradients, vorticity shedding, and localized small length scale flow features, dynamic adaptive mesh refinement (given by hierarchical nested Cartesian grid patches) is employed. That spatial dynamic adaptation is used in conjunction with a variable time step, second-order, linearly implicit time integration scheme. The capabilities of the present numerical method are demonstrated by numerical verification and simulation of two classical diffusion flame examples from literature.

Simulation of the inelastic deformation of porous reservoirs under cyclic loading relevant for underground hydrogen storage

Subsurface geological formations can be utilized to safely store large-scale (TWh) renewable energy in the form of green gases such as hydrogen. Successful implementation of this technology involves estimating feasible storage sites, including rigorous mechanical safety analyses. Geological formations are often highly heterogeneous and entail complex nonlinear inelastic rock deformation physics when utilized for cyclic energy storage. In this work, we present a novel scalable computational framework to analyse the impact of nonlinear deformation of porous reservoirs under cyclic loading. The proposed methodology includes three different time-dependent nonlinear constitutive models to appropriately describe the behavior of sandstone, shale rock and salt rock. These constitutive models are studied and benchmarked against both numerical and experimental results in the literature. An implicit time-integration scheme is developed to preserve the stability of the simulation. In order to ensure its scalability, the numerical strategy adopts a multiscale finite element formulation, in which coarse scale systems with locally-computed basis functions are constructed and solved. Further, the effect of heterogeneity on the results and estimation of deformation is analyzed. Lastly, the Bergermeer test case-an active Dutch natural gas storage field-is studied to investigate the influence of inelastic deformation on the uplift caused by cyclic injection and production of gas. The present study shows acceptable subsidence predictions in this field-scale test, once the properties of the finite element representative elementary volumes are tuned with the experimental data.

Post-buckling analysis of shell-like structures using an implicit dynamic time integration scheme

PurposeThis paper aims to investigate post-buckling responses of shell-like structures using an implicit conservative-decaying time integration dynamic scheme.Design/methodology/approachIn this work, the authors have proposed the use of a four-node quadrilateral flat shell finite element with drilling rotational degree of freedom within the framework of an updated Lagrangian formulation mutually with an implicit conservative-dissipative time integration dynamic scheme.FindingsSeveral numerical simulations were considered to evaluate the accuracy, robustness, stability and the capacity of the considered time integration scheme to dissipate numerical noise in the presence of high frequencies. The obtained results illustrate a very satisfying performance of the implicit conservative-dissipative direct time integration scheme conjointly with the quadrilateral flat shell finite element with drilling rotation.Originality/valueThe authors have investigated the potential of the implicit dynamic scheme to deal with unstable branches after limit points in the non-linear post-buckling response of shell structures with no need for structural damping. The capability of the studied algorithm to study buckling and post-buckling behaviour of thin shell structures is illustrated through several numerical examples.

Imposing displacements in implicit direct time integration & a patch test

We consider the transient analysis of finite element models when an implicit direct time integration scheme is used and focus on the accuracy of the response solution when time varying nodal displacements are imposed. The solutions of such problems are abundantly sought in engineering practice. We analyze the performance of the widely used Trapezoidal Rule, a specific case of Newmark time integration, and the performance of the Bathe method. The theoretical results show that the Trapezoidal Rule is unstable and this instability becomes more pronounced as the time step size for solution is decreased, whereas the Bathe method is stable. We illustrate these theoretical findings in the transient analyses of some finite element models. These investigations lead us to define a “patch test” useful for the general evaluation of direct time integration schemes.